Tracking energy fluctuations from fragment partitions in the Lattice Gas model 
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Partial energy fluctuations are known tools to reconstruct microcanonical heat capacities. For ex- 
perimental applications, approximations have been developed to infer fluctuations at freeze out from 
the observed fragment partitions. The accuracy of this procedure as well as the underlying inde- 
pendent fragment approximation is under debate already at the level of equilibrated systems. Using 
a well controlled computer experiment, the Lattice Gas model, we critically discuss the thermody- 
namic conditions under which fragment partitions can be used to reconstruct the thermodynamics 
of an equilibrated system. 



I. INTRODUCTION 

The observation of abnormal fluctuations in nuclear 
multifragmentation and its possible connection to a neg- 
ative heat capacity has raised much interest and dis- 
cussions in the last years0ISSISIE0- 

For a microcanonical ensemble, it has been proposed in 
ref. || that the kinetic energy fluctuation cr\/T 2 can be 
used to reconstruct the heat capacity even in the context 
of phase transitions and for small systems. The accuracy 
of the fluctuation expression has been successfully tested 
on numerical experiments on the liquid gas phase transi- 
tion, using the microcanonical Lattice Gas model[|| and 
molecular d yna mics simulations with a Lennard-Jones 
potential |llTol|. 

From the experimental point of view, it has been pro- 
posed to use the clusters asymptotically detected in nu- 
clear multifragmentation reactions to backtrace the fluc- 
tuations of the total energy partitioning at freeze-out 
p|. The robustness of the experimental procedure has 
been tested in ref. ,11]. In particular, statistical models 
have been used to generate events, then analyzed using 
the experimental procedure and a good reproduction of 
the model heat capacities has been reported Us- 
ing molecular dynamics at equilibrium and recognizing 
fragments through the Hill algorithm[l2| the authors of 
ref. [j| have come to opposite conclusions criticizing the 
independent fragment hypothesis on which is based the 
experimental method as well as the zero temperature ap- 
proximation for the fragment binding used to evaluate 
the fragment internal interaction energy. The model of 
ref. has also been studied in ref. 07 where the dynam- 
ics of the expansion was explicitly included. The result 
was that the fluctuations of dense systems were strongly 
modified by the dynamics and that only the latest phase 
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of the expansion, after a freeze-out corresponding to di- 
lute configurations, can be accessed from the observation 
of asymptotic partitions. 

More generally, the question of the thermodynamic 
conditions under which the energy partitioning of a small 
fragmenting system can be retraced from the measured 
fragment sizes and kinetic energies, raises important 
questions for the whole field of nuclear thermodynamics. 
Indeed the independent fragment hypothesis is not only 
needed to reconstruct partial energy fluctuations but 
is also necessary for any other quantitative estimation of 
the nuclear phase diagram^]. Moreover only if the frag- 
menting source can be approximated by an ensemble of 
non- (or weakly-) interacting nuclear clusters, the statis- 
tical models [13 that have successfully reproduced heavy 
ion data since two decades can be theoretically justified. 

To contribute to this debate, we want to address the 
problem of the independent fragment hypothesis and of 
the fragment energetics at equilibrium in the framework 
of a well controlled exact numerical model, the Lattice 
Gas model(3- We will show in this article that, in the 
Lattice Gas model, the system heat capacity can be 
well estimated from fragment sizes at all energies and 
for all pressures p/p c ^$ 1/3 almost independent of the 
parametrization adopted for the fragment energies. For 
higher pressures this approximation tends to break down, 
but the estimation of fluctuation stays at the 30% accu- 
racy level even in the supercritical regime. 



II. LATTICE GAS RESULTS 

It has been proposed in to use the partitioning of 
a fixed total energy (E tot ) into kinetic (K) and configu- 
rational (V) energies, E to t = K + V, in order to look for 
entropy curvature anomalies. Indeed, for classical sys- 
tems with momentum independent interactions, because 
of the microstates equiprobability the microcanonical dis- 
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tribution of K at a fixed energy E tot reads 



liquid drops contributions 



P E (K) = cxp (S K (K) + S v (V) - S tot (E tot )) 

where Sk, Sy, S to t are the kinetic, potential, and total 
entropies. If this distribution is normal, a gaussian ap- 
proximation can be performed leading to analytic expres- 
sions relating temperatures and heat capacities to the 
observed moments of the distribution. When only the 
leading order is kept, we recover the simplest expression 
used in ref. 0] to relate the total microcanonical heat 
capacity C to the kinetic one Ck and the kinetic energy 
fluctuation ok scaled by the system temperature T 
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Even if in practice it seems that Eq. is often accurate 
enough, the validity and the accuracy of this approxi- 
mate expression should always be checked by controlling 
the actual distribution and when the distribution is nor- 
mal by evaluating the corrections. In our work this has 
been done both in theory^ and experiments [TO] . The 
accuracy of the simple expression Q to infer the heat 
capacity even in the vicinity of a phase transition, or 
worse of a critical point is also verified in refs. 0, . 

The application to experimental data of the idea pro- 
posed in [H requires the development of different tools. 
Let us first briefly recall the problem of the potential 
energy determination for an ensemble of fragments. If 
we look at a system of A interacting particles as a 
system of Mf clusters (including monomers), the po- 
tential energy V = 



EjL'i Vf + Ef< 9 V fa where V f = £ 



can be written as V — 



i<jef v n 



and V fg 
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are the intrafragment and interfragmcnt 
components respectively. In the experimental analysis of 
ref.[l| the only interaction considered among the different 
fragments is the Coulomb force because of its long range 
nature. The interfragment nuclear force is thus neglected 
following the argument that an important nuclear inter- 
action is incompatible with the freeze-out concept. Con- 
cerning the evaluation of the intrafragment potential en- 
ergy Vf = Ei<je/ v ij> this latter is approximated in the 
experimental analysis^ by the tabulated ground state 
nuclear energies. Both these approximations are used in 
man y o ther thermodynamic studies of multifragmenta- 
tion |l4llTs| and in all macroscopic statistical mo dels [TH- 



A. Independent Liquid drop approximation 



qld = J2 B f = Y1 avA f + asA ) 



2/3 



(2) 



where Af is the size of cluster / and Mf is the total 
number of clusters. The parameters ay, as can then 
be fitted to reproduce the cluster binding at low but fi- 
nite temperature to avoid the peculiarities of the ground 
states of such classical systems. 

In this article we present a study based on exact nu- 
merical experiments of the 3D Lattice gas model 9j . The 
ground states of such a model are cubes and, as discussed 
above, to estimate the internal interaction energy of the 
liquid fragments avoiding the extra binding of those pe- 
culiar cubic configurations we shall use low temperature 
simulations. Using canonical simulations with temper- 
atures around 1/3 of the critical temperature (or mi- 
crocanonical ones with energies around — 1.5e) leads to 
ay = — 2.86e, as = 2.73e, where e is the lattice coupling 
(see discussion of Fig 1). 

For comparison the choice ay — — 3e, 05 = 3e cor- 
responds to large cubic clusters, while ay = — 3.06e, 
as = 3.35e leads to a good description of the zero tem- 
perature clusters in the size range 2 < Af < 60. 

In the following calculations the coefficients ay = 
— 2.86e, as = 2.73e will be kept constant, and we will 
come back to the influence of the parameters values in 
the last section. 
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In order to check the quality of these approximations, 
we shall use an exact numerical experiment. This re- 
stricts studies to classical systems. At variance with nu- 
clear systems which are liquid in their ground states, the 
ground states of classical models are solid and so present 
an extra binding. To avoid this difficulty, we may approx- 
imate the interaction energy V as a sum of independent 



FIG. 1: Correlation between the exact interaction energy V 
and its approximation from fragment sizes Qld (see text) at 
an average density p/ po = 0.0135 and a temperature T/T c = 
0.29 (upper part), and T/T c = 0.57 (lower part). The liquid 
drop parameters for the fragment binding energy are fixed as 
a v = -2.86e, a s = 2.73e. 
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coefficient is c = 0.77 ) shows that indeed in such a low 
density configuration the interaction energy can be cal- 
culated within a liquid drop approximation. 

Even in the independent fragment picture, the liquid 
drop coefficients are expected to change with temper- 
ature reflecting the internal excitation of the produced 
clusters. This can be seen in the lower part of figure 
^ which shows a calculation for the same very diluted 
system at a temperature well above the transition tem- 
perature. It is clear that the quality of the correlation 
does not significantly change with the temperature (here 
the Qiu/y-correlation coefficient is c = 0.78), but the 
use of the low temperature mass formula leads to a sys- 
tematic Ri 20% overestimation of the fragment binding 
that can affect the thermodynamic analyses. 



FIG. 2: Reconstruction of the average kinetic energy (solid 
lines) as a function of total energy from fragment sizes in 
the Lattice Gas model at different pressures. Dashed lines: 
liquid drop parameters from table 1. Symbols: liquid drop 
parameters fixed from the low temperature low density phase 
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FIG. 3: Reconstruction of the normalized kinetic energy fluc- 
tuation (solid lines) as a function of total energy from frag- 
ment sizes in the Lattice Gas model at different pressures. 
Dashed lines: liquid drop parameters from table 1. Symbols: 
liquid drop parameters fixed from the low temperature low den- 
sity phase as in figure 1. 



We have first tested the accuracy of the approximation 
for the total interaction energy V in a low density 
case (108 particles in a cubic box of volume L 3 = 8000) 
where the interfragment energies can be safely neglected. 
The upper part of figure ^ shows the correlation between 
the exact interaction energy and its independent frag- 
ment approximation eq. Q in a canonical simulation at a 
temperature T/T c = 0.29. The corresponding average to- 
tal energy is — 1.74e, well below the backbending region. 
The good linearity of the plot (the QLo/I^-correlation 



B. The choice of the statistical ensemble 

This point should be further explored looking at di- 
rect effects of the considered approximation on thermo- 
dynamical quantities, i.e. on ensemble averaged observ- 
ables such as averages and variances. Since the differ- 
ent statistical ensembles are not equivalent in finite sys- 
tems, the statistical ensemble has to be specified. In 
the following we have chosen to perform calculations 
in the microcanonical "isobar" ensemble characterized 
by the two state variables (E, A), the total energy and 
the Lagrange multiplier imposing the average volume V, 
respectively 9] . Statistical averages are calculated as 



J2 (n) A^exp (-j3E^ - WW) 5 (E^ - E) 
ex P (-^ (n) - AV(™)) S (£W - E) 

(3) 

where A is a generic observable (A = K allows to com- 
pute the microcanonical temperature, A = K 2 provides 
the heat capacity), the sum runs over the system mi- 
crostates, and the average volume is defined through the 
one body observable 



i—1 



(4) 



where rii = 0, 1 is the occupation of the i — th lattice 
site. 

Different reasons motivate the choice of this ensemble. 
First, this is the ensemble in which the liquid gas phase 
transition is associated to a negative heat capacity up 
to the critical point [19(. Moreover in the actual analysis 
of heavy ion experiments data are sorted in excitation 
energy bins, i.e. approximate realizations of microcanon- 
ical ensembles. On the other side the system extension is 
only imposed by the freeze out requirement without any 
boundary condition. We have recently shown that such 
a physical situation is accounted in a thermostatistical 
coherent way only if the system size is imposed through 
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a lagrange parameter A, with P = X/ (3 having the dimen- 
sion (and the physical role) of a constraining pressure po). 
Finally it is interesting to notice that if A = this "iso- 
bar" ensemble is equivalent to the isochore ensemble used 
in section lTl Al if the volume is large enough such that the 
boundary conditions become irrelevant. 



C. Averages and fluctuations 

Figure [21 shows the estimated average kinetic energy 
E — J2f Bf a,s & function of the total energy with differ- 
ent constraints on the system volume :9], corresponding 
to different pressures. The shape of the exact caloric 
curve (solid line) is nicely reproduced by the indepen- 
dent "cold" liquid drop approximation (open symbols) 
for all volume constraints. Using the temperature back- 
bending to define the phase transition, we see that the 
coexistence zone as well as the critical point can be well 
estimated from the unique knowledge of the fragment 
partitions. However, the actual value of the tempera- 
ture shows a systematic shift. The fact that this shift is 
also present in the pure high temperature low density gas 
phase (figure suggests that it may be due to a temper- 
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FIG. 4: Upper panel: partial energies of an ensemble of events 
at a pressure X/X c = .33 and energy .38e < E to t/N < A2e 
in the middle of the coexistence region. From bottom to top: 
-interaction energy per particle; -its independent fragment ap- 
proximation (LD parameters as in figure 1) (shifted by 0.4); 
-total energy per particle. The horizontal lines give the aver- 
age values of the corresponding energies. Lower panel: same 
as the upper panel on an extended scale. 

The effect on partial energy fluctuations is shown in 
figure |21 for the same thermodynamic conditions as in 
figure El We can see that the fluctuations tend to be 
overestimated, but since the bias on the average value 
goes in the same direction, the normalized fluctuations 



are still reasonably reproduced almost up to the critical 
point. 
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FIG. 5: Qld /V -correlation coefficient (upper part) and per- 
centage error \oq — o\\la\ on the partial energy fluctuation 



from figure (lower part). The liquid drop parameters are 
fixed as av ~ — 2.86e, as = 2.73e. Black points:X/X c = .04, 
squares:\/ \ c = .33, triangles :X / X c = 1, stars: X/X c — 2.46. 

Even if the quantitative study of a statistical ensemble 
should be based on ensemble averages such as the one 
presented in the previous figures, it is instructive to look 
at a set of events corresponding to one of the above cases 
as shown in Fig. 0] This figure illustrates the fact that 
the actual interaction energy and its independent liquid- 
drop approximation have similar properties, averages and 
fluctuations. Moreover, the extended scale part shows 
that the two quantities are correlated. 



D. Errors and correlations 

In order to quantify the quality of the liquid drop ap- 
proximation let us look at ensemble averaged quantities 
which are indeed the information used to infer thermody- 
namical properties. One of the most important observ- 
ables is the kinetic energy fluctuation. The correspond- 
ing percentage error \oq — oWja 2 K is shown in the lower 
part of figure One observes that the error increases 
with the pressure but remains low for all energies up to 
the critical pressure. This error is below 10% up to 1/3 of 
the critical pressure and below 20% at the critical pres- 
sure. Moreover, when used to compute the heat capacity 
this error is partly compensated by the similar error on 
the kinetic energy used to deduce the temperature. This 
explains the good accuracy observed for the liquid drop 
approximation up to rather high temperatures and pres- 
sures. 

This good accuracy of the ensemble averaged quanti- 
ties does not imply nor require a similar accuracy on an 
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event by event basis. In fact one expects stronger devia- 
tions when looking at a single event. One way to study 
this event by event accuracy of the liquid drop approxi- 
mation is to study the correlation coefficient c. However, 
this is only a part of the discussion since the correla- 
tion coefficient is not sensitive to the magnitude of the 
fluctuation and only characterizes the link between the 
independent liquid drop approximation and the actual in- 
ternal energy. In the considered case, this implies that by 
construction the correlation coefficient c is independent 
of the value of the liquid drop parameters. The global 
trend of the correlation coefficient with pressure and en- 
ergy is presented in the upper part of figure [5] We can 
see that the best correlation is systematically observed 
around the fluctuation peak, where the size distribution 
is the broadest. At subcritical pressures, c always ex- 
ceeds 0.4, while the correlation decreases to around 0.2 
in the supercritical regime. 

It should be stressed that the idea to use, in equation 
(1), the first moments (averages and fluctuations) of the 
event distribution was indeed to overcome the difficulty 
to get accurate information on an event by event basis. 
Though the best correlation obviously leads to the small- 
est error, we can see that the relation between c and Act 
is not trivial. In particular a value of c exceeding 0.4 
comes out to be enough to keep the error below the 15% 
level, independent of the energy. 

It is also interesting to remark that the correlation co- 
efficients shown in figure [S] are systematically lower than 
the ones shown in the calculations of figure 2] This is 
another illustration of the non equivalence of statistical 
ensembles: in the canonical case the distributions are 
wider than in the microcanonical one allowing a better 
correlation. 



E. Beyond the fixed cluster energy approximation 

The failure at increasing pressure shown by fig|3|is in- 
teresting. One may wonder whether the lack of reproduc- 
tion is due to the breaking down of the independent frag- 
ment approximation |5j in dense media, or whether the 
configurational contribution to the internal fragment ex- 
citation energy has to be taken into account by a proper 
redefinition of the mass formula. To answer to this ques- 
tion we have allowed a free variation of the liquid drop 
parameters according to table 1. The result, shown by 
the dashed lines in figs l2!3l is that both the caloric curve 
and the fluctuations can be very precisely reproduced by 
the independent fragment approximation in a wide range 
of temperatures and pressures if the volume and surface 
coefficients are allowed to decrease with increasing ex- 
citation. Only at very high pressures and temperature, 
above the critical point, the approximation appears to 
clearly break down: in this dense configuration the ob- 
jects identified as "fragments" by the cluster recognition 
algorithm have certainly little to share with physical iso- 
lated liquid drops. Indeed in this density regime the Q 



value eq. is very poorly correlated with the interaction 
energy (see figure EJ). 
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TABLE I: Volume ay and surface as effective liquid drop pa- 
rameters (see text) allowing to reproduce the average canoni- 
cal configurational energy of 216 Lattice Gas particles at the 
temperature corresponding to the maximal energy fluctuations, 
for different pressures normalized to the critical pressure. The 
average volume occupied by the system divided by the critical 
volume is also given. 



III. DISCUSSION AND POSSIBLE 
EXPERIMENTAL IMPLICATIONS 

Let us now turn to the possible implications of these 
findings to heavy ion collisions experiments, and partic- 
ularly to the determination of the heat capacity from the 
fluctuations of asymptotic detected partitions. 

In the application of eq^to the analysis of multifrag- 
mentation data several difficulties arise that have to be 
considered. 

First of all in the experimental case the only infor- 
mation available on the system is given by its clus- 
ter properties (sizes and kinetic energies). This means 
that the total energy of the system is not evaluated 
from its microscopic constituents, but is also estimated 
from the measured kinetic energies and fragment sizes as 
E tot = Qld + K. 

Moreover both Qld and K are time dependent vari- 
ables, and in the data analysis the quantities at freeze out 
are extrapolated from the asymptotic ones solely correct- 
ing for secondary evaporation and Coulomb repulsion. In 
particular the configurational energy functional associ- 
ated to each fragment at freeze out is assumed to be the 
one associated to T = and is not pressure dependent. 
As we have seen in the previous section, this assumption 
is not valid for the Lattice Gas model: the presence of a 
cubic ground state in this model leads to a mass formula 
that is not adequate to reproduce the average cluster en- 
ergy in the liquid phase. The difference in energy is small 
compared to the liquid- vapour latent heat, however we 
have seen that it can have important effects in the cal- 
culation of the observables. This difficulty in principle 
should not arise in nuclear physics where clusters are al- 
ready liquid at T=0, however one may ask how much a 
possible modification with temperature and pressure of 
the fragments energetics would influence the experimen- 
tal analysis. 

Finally the assumption is explicitly made that at the 
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freeze out time fragment partitions essentially reflect 
thermal equilibrium. 
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FIG. 6: Same as figure^ above, but including the effect of the 
calorimetric estimation of total energy (see text). Squares: 
zero temperature liquid drop parameters a v — —3e,a s = 3e. 
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FIG. 7: Same as figure^ above, but including the effect of the 
calorimetric estimation of total energy (see text). Squares: 
zero temperature liquid drop parameters a v — — 3e,a s = 3e. 



A. The thermodynamics of a gas of clusters 

Let us leave the difficult problem of equilibrium and 
time dependence to the next subsection and start ad- 
dressing the question of the quality of reconstruction of 



the thermodynamics of freeze-out, using accessible infor- 
mation (fragment sizes and total kinetic energies). In 
this matter, equilibrium studies provide a valuable test- 
ing ground, if the studied thermodynamic conditions cor- 
respond to the freeze-out configurations of the model. 

To quantify the uncertainty of the heat capacity recon- 
struction from the fragment information only, we have 
introduced reconstructed microcanonical statistical en- 
sembles defined by the total energy constraint E ca i — 
K + Qld at different pressures and with different pre- 
scriptions for the liquid drop parameters. Figures|S]and[7| 
show the resulting first and second moment of the K dis- 
tribution in bins of E ca i . We can see that the systematic 
bias observed in fig|21 disappears if the same prescription 
for the interaction energy is used both for the fluctuation 
and for the total energy of the system. Up to the critical 
point, the precise parametrization of the liquid drop en- 
ergy does not change the results dramatically. In partic- 
ular the zero temperature cubic solid cluster prescription 
a v = — 3e,a s = 3e induces a spurious shift towards higher 
energies, but does not induce extra fluctuations. 

These results may be qualitatively understood as fol- 
lows. The error in the estimation of V through Qld in 
each event of a set at a given total physical energy in- 
duces an overall bias in the evaluation of < E to t >, but 
also and more important, a spread in the total energy 
estimation. Both these effects are especially important 
at high pressures, where the correlation coefficient be- 
tween Qld and V is low (see figE}- The systematic bias 
in the estimation of the average total energy depends di- 
rectly on the value chosen for the liquid drop parameters, 
and leads to the shift towards higher energy observed in 
figs lul7l The width of the E ca i distribution for each value 
of E to t is at the origin of the overestimation of the par- 
tial energy fluctuations observed in figlHl When events 
are analyzed in bins of total estimated energy, these spu- 
rious energy fluctuations do not contribute any more to 
the width of the K distribution. The remaining effects 
are due to the mixing in a given E ca i bin of events coming 
from different physical E to t values. Since this mixing is 
approximately symmetric in energy, the average < Ek > 
appears to be not much affected and the main effect of 
mixing is to flatten out the fluctuation curve. This is 
very similar to the results already reported for the anal- 
ysis of events produced in macroscopic statistical models 
(SMMpI and SIMON^U). 

From figs and [7| we can also see that, when calorime- 
try is taken into account, it is not possible any more to 
tune the liquid drop parameters at high pressure such as 
to reproduce at the same time the correct average energy 
and fluctuation. This confirms again that the indepen- 
dent fragment hypothesis breaks down in the supercriti- 
cal regime. 

The error on the fluctuation estimation is reported in 
figEl f° r all t ne considered pressures and energies, using 
the fixed low temperature liquid drop parameters ay = 
— 2.86e, as = 2.73e. The fact of taking into account the 
calorimetric shift does not increase this systematic error 
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FIG. 8: Percentage error |<tq — ct^-I/ct^ on i/ie partial en- 
ergy fluctuation from figure\T\ The liquid drop parameters are 
fixed as av = — 2.86e, as = 2.73e. Black points:X/X c = .04, 
squares:\/ X c = .33, triangles :X/X C = 1, stars: X/X c = 2.46. 

above the 20% level at subcritical pressure. 

The results of figs l6l7l suggest that if the freeze-out 
density is low enough, then the thermodynamic prop- 
erties of the system can be deduced from the fragment 
partitions using the independent fragment model with 
a fragment energy approximated by its low temperature 
binding. In the Lattice-gas model the validity condition 
of this simple approximation happens to be not so restric- 
tive since it corresponds to an average volume V > 1.5V C . 
If the freeze out density is higher, in the lattice-gas model 
the independent fragment approximation tends to break 
down, and both the temperature and the fluctuation are 
systematically underestimated. 



B. The influence of dynamical and quantum effects 

When dealing with reactions and trying to extract 
thermodynamic information, the first question is the rel- 
evance of the freeze-out and equilibrium concepts at spe- 
cific stages of the collision. From the theoretical point of 
view, this complex question, which is strongly debated 
since more then a decade, can only be addressed by com- 
plete simulations of the reaction and critical analysis of 
the resulting time dependence. Such studies require to go 
beyond the equilibrium models used in the present article 
as well as in ref. p|. In a recent dynamical calculation 
with the Lennard-Jones hamiltonian, A.Chernomoretz et 
al.[To| have addressed this question by looking at the time 
evolution of different observables for a system initially 
thermalized in a confined box and subsequently freely 
evolving in the vacuum. Almost independent of the ini- 
tial energy and density, the average configurational en- 



ergy and the associated fluctuations turn out to freeze 
when the density of the system is of the order of p « 0.03 
in Lennard-Jones units, a region of the phase diagram 
well below the critical point |2j. This study, thus, im- 
plies that only such diluted stages of the reaction can 
be reconstructed from the asymptotic ("experimental") 
information. However, it is important to stress that the 
density value at freeze out is a model dependent quan- 
tity and in particular it is correlated with the range of 
the force. In the study of ref. [lCf . the configurations 
at freeze-out are also characterized by an interfragment 
energy sufficiently small S/< S V/ S (i) ~ f° r fragment 
partitions to be essentially frozen. Conversely to the ab- 
solute value of the density, this condition might be a more 
robust definition of freeze out. It is also a first argument 
in favor of an independent fragment approximation. 

Another consequence of the time dependence of a col- 
lisional process is that the freeze out configuration may 
not be sufficiently close to a thermodynamic equilibrium 
to be described with statistical toolsy, H3| • The statis- 
tical nature of freeze-out is indeed the key point in or- 
der to interpret the reconstructed fluctuation in terms of 
heat capacity. To estimate the distortions due to the out 
of equilibrium component one needs to know how much 
the averages and fluctuations deviate from the equilib- 
rium values. Different verifications have been performed 
to estabilish this point [Tl| in the experimental analyses 
of collisional data. In particular, the observation of the 
same behavior with different entrance channels is an in- 
teresting argument in favor of being close to an equilib- 
rium. From a theoretical point of view, this fundamental 
open question can only be addressed through dynami- 
cal approaches and cannot be answered from statistical 
calculations as the one presented here. 

Finally it is important to stress that the results of clas- 
sical models like the Lattice Gas (or the Lennard Jones 
analyzed in refs. cannot be quantitatively ap- 

plied to nuclear data because of the complete lack of 
quantum effects. The ground state properties are not 
the only point where quantum effects are expected to 
play an important role. The actual value of the density 
at freeze out of a nuclear system can be very different 
from the one estimated from classical calculations [lOj, 
and the same is true for the limiting density that can 
be described through the independent fragment approxi- 
mation derived in this work. The consistent inclusion of 
quantum effects in the statistical analyses of nuclear data 
is an ambitious program that is only in its infancy ^ij 
and constitutes one of the greatest theoretical challenges 
of heavy ion collisions in the next decades. 



IV. CONCLUSION 

To summarize, in this paper we have reconstructed 
configurational energy fluctuations from the fragment 
partitions of a Lattice Gas model. As far as the first and 
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second moments of collective observables are concerned, 
dilute systems in thermodynamic equilibrium can be ac- 
curately approximated by an ensemble of independent 
fragments. A unique liquid drop parametrization for the 
fragment binding energies is able to describe the ther- 
modynamics of the system independent of the deposited 
energy or temperature. A modification of the order of 
10% of the liquid drop parameters does not modify the 
temperature and fluctuation in a sizeable way, if the same 
parametrization is consistently employed for the determi- 
nation both of the fluctuation and of the total energy. For 



pressures of the order or above the critical point, the in- 
dependent fragment approximation tends to break down 
leading to a systematic underestimation of temperatures 
and fluctuations at the 30% level. 
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